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ABSTRACT 

Context. We report the detection of oxygen-rich circumstellar envelopes in stars of the nearby (700 kpc) starburst galaxy IC 10. The 
star formation history and the chemical environment of this galaxy makes it an ideal target to observe dust production by high-mass 
stars in a low-metallicity environment. 

Aims. The goal of this study is to identify oxygen-rich stars in IC 10 and to constrain their nature between asymptotic giant branch 
stars (AGBs), red supergiants (RSGs), and other infrared bright sources. We examine the mass-loss rate of the stars and compare to 
results obtained for the Magellanic Clouds. Our objectives are to (1) assess whether RSGs can be significant dust producers in IC 10, 
and (2), solve the discrepancy between the star formation history of IC 10 and the relatively low number of RSGs detected in the 
optical. 

Methods. We search for silicate dust in emission by using the spectral map observed with the Infrared Spectrograph on board the 
Spitzer Space Telescope. The optical (UBVRI) and infrared (JHK, Spitzer/IRAC and Spitzer/MIPS) photometry is used to assert the 
membership of the stars to IC 10 and disentangle between AGBs and RSGs. Radiative models are used to infer mass-loss rates and 
stellar luminosities. 

Results. The luminosity and colors of at least 9 silicate emission sources are consistent with stars within IC 10. Furthermore, the 
photometry of 2 of these sources is consistent with RSGs. We derive dust mass-loss rates similar to the values found in the Magellanic 
Clouds. Accounting for the sample completeness, RSGs are not important contributors to the dust mass budget in IC 10. 

Key words. Stars: mass-loss - Galaxies: individual :IC 10 - Infrared: stars 



1. Introduction 

Principal contributors to the dust in the interstellar medium 
(ISM) are thought to be low-mass asymptotic giant branch 
(AGB) stars, supernovae (SNe), red supergiants (RSGs), and 
late-type WC Wolf-Rayet stars (e.g., Gehrz 1989). AGB stars 
are expected to contribute to significant dust production in the 
most metal-poor sources, because their evolution timescales 
are shorter in such environments. For metallicities as low as 
1/2OOZ , stars might take as short as lOOMyr to evolve from 
the zero-age main sequence to the AGB (Ventura et al. 2002; 
Herwig 2004). Boyer et al. (2012) found that AGBs produce 
most of the dust from cool evolved stars in the metal-poor 
Small Magellanic Cloud (SMC, 0.2 Z Q ), but other dust produc- 
tion sources such as growth from existing grains are necessary to 
explain the total dust mass budget. Alternatively, it is largely de- 
bated whether cosmic dust abundance can be reconciled with SN 
dust (e.g., Li et al. 2008; Maiolino et al. 2004), as dust mass pro- 
duced in SNe is largely uncertain (10~ 3 M up to 0.5 M Q ; e.g., 
Stanimirovic et al. 2005; Matsuura et al. 2011). With these con- 



siderations in mind, high-mass stars could be significant con- 
tributors to the ISM dust. In metal-poor galaxies, where late- 
type WCs are scarce, RSGs might be the dominant dust source 
(Massey et al. 2005; but see Boyer et al. 2012). 

The detection of infrared excess due to mass loss in the cir- 
cumstellar envelopes of evolved stars within Local Group galax- 
ies is mostly limited to low-mass carbon-rich stars (~ 1 M©; e.g., 
Sloan et al. 2012 and references therein). This is due to the star 
formation history of our nearest neighbor galaxies, with the lack 
of a recent starburst episode that would result in pronounced 
populations of - more massive - oxygen-rich stars. From an 
observational point of view, probing O-rich dust production in 
Local Group galaxies is challenging as it requires observing the 
silicate emission bands in the mid-infrared range. The observa- 
tion of such stars has been therefore limited only to the Milky 
Way and to the Magellanic Clouds with the space telescopes ISO 
(e.g., Trams et al. 1999a; 1999b) and Spitzer (e.g., Buchanan et 
al. 2009; Groenewegen et al. 2009; van Loon et al. 2010; Woods 
et al. 2011; Boyer et al. 2011). The great sensitivity of Spitzer 
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along with matured data analysis techniques now makes it possi- 
ble to observe O-rich dust spectra beyond the Magellanic Clouds 
and in starburst galaxies. 

The nearby dwarf starburst IC 10 is an ideal target to detect 
O-rich stars in a metal-poor environment, such as what has been 
accomplished for the Magellanic Clouds. IC 10 was discovered 
by May all (1935) and Hubble (1936); it is the nearest starburst 
galaxy known (-700 kpc; Kim et al. 2009; Kennicutt et al. 1998; 
Borissova et al. 2000; Hunter 2001). Its size and mass are com- 
parable to the SMC while its metallicity of 12 + log(0/H)^ 8.26 
(1/2.7 Z assuming the solar abundance from Asplund et al. 
2009) lies between that of the Small and Large Magellanic Cloud 
(Garnett 1990; Lequeux et al. 1979; Richer et al. 2001 ; Skillman, 
Kennicut & Hodge 1989). The starburst nature of IC 10 was first 
revealed by the discovery of a large number WR stars by Massey 
et al. (1992). IC 10 has experienced several episodes of exten- 
sive star formation, with the most recent ones a few 10s to 100s 
of Myr ago (Vacca et al. 2007). The presence of a widespread 
population of WR stars (Massey et al. 1992; Massey & Homes 
2002; Royer et al. 2001) suggests that IC 10 starburst is also 
widespread and that RSGs are to be expected throughout the 
galaxy. At low metallicities, the number of RSGs should even 
dominate over the number of WR stars, as one expects from 
evolutionary theory (Maeder et al. 1980) and as demonstrated 
observationally in the Local Group by Massey (2002, 2003). The 
number of spectroscopically confirmed IC 10 WRs is 24, and the 
actual number is believed to be many more (Massey & Holmes 
2002). From the relatively low metallicity of IC 10, one would 
thus expect the population of RSGs to be about 50 - 100 strong 
(see Figure 12 of Massey 2003), not inconsistent with the color- 
magnitude diagram (CMD) of this galaxy. 

After presenting the observations in Section[2] we derive 
a preliminary silicate emission map of IC 10 from which we 
identify several O-rich candidates (Section[3). We then cross- 
correlate our sample with optical and infrared catalogs, remov- 
ing foreground stars in the process (Section]?]). The final mid- 
infrared spectra (Section]?]) are used to derive mass-loss rates 
and discuss the stellar chemistry of the sources (Section|6]). 



2. Observations 

IC 10 was observed with the Infrared Spectrograph (IRS; Houck 
et al. 2004) onboard the Spitzer Space Telescope (Werner et al. 
2004) on 2008 September 13 as part of GTO program 50318. 
Observations consisted in a sparse spectral map with the Short- 
Low (SL) module, providing the wavelength coverage between 
« 5 - 14.5 jim with a spectral resolution power between 60 and 
130. The map is made of 58 perpendicular steps and 8 parallel 
steps, with 2 cycles of 14 seconds per exposure. The galaxy was 
not fully sampled spatially because of time constraints. A gap 
was deliberately introduced between every perpendicular scan- 
ning position, with the gap size precisely equal to the width of 
the Short-Low aperture 3.7"). The full width at half maxi- 
mum (FWHM) of the point spread function is on the order of 
the aperture width (« 3.5" at lOyum), allowing the detection 
of light from outside the slit despite the gaps. The flux is cali- 
brated by performing optimal extraction of the point-like sources 
and accounting for the slit throughput (Section]?]). In order to 
solve the incomplete spatial sampling for the preliminary analy- 
sis (Section[3|, gaps were interpolated using a cubic spline over 
2 pixels on each side. We estimate the effective spatial resolution 
of the map to be somewhat larger than the resolution at 14.0 yum, 
i.e., FWHM> 1.8 px = 3.7". 



A preliminary automatic cleaning of each exposure was per- 
formed using irscleanQ The data were then imported and an- 
alyzed with CUBISM (version 1.7; Smith et al. 2007). A second 
manual cleaning step was performed using the backtracking tool 
provided by CUBISM. Several exposures at the edges of the map 
were chosen to remove the background emission, mostly arising 
from the Milky Way. Images corresponding to relevant wave- 
length ranges for building the silicate strength map (Section[3]) 
were extracted with CUBISM. The map of the integrated flux in 
the mid-IR range is shown in Figure]!] 

Near-infrared photometry of stars toward IC 10 is taken 
from the 2MASS point-like source catalog (Skrutskie et al. 
2006) while optical photometry is taken from the Local Group 
Galaxies Survey (LGGS; Massey et al. 2007). In addition, we 
measured the mid-IR photometry using IRAC (Fazio et al. 2004) 
and MIPS (Rieke et al. 2004) onboard Spitzer. IC 10 was ob- 
served with IRAC on 2004 July 23 with the 4 channels, centered 



1 Version 1.9; http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 




40.0 30.0 20.0 10.0 0:20:00.0 19:50.0 



Right Ascension (J2000) 

Fig. 1. Spitzer/IRS map. The top panel shows the distribution 
of the flux integrated along 7.4 - 14.5 //m. The effective spatial 
resolution of the map (> 3.2") is approximately equal to the 
width of the SL aperture. The silicate strength map is presented 
in the bottom panel (Section|3]). 
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respectively at 3.6 /mi, 4.5 fim, 5.8 yum, and 8.0/mi (AORKEY 
4424960), and with MIPS on 2004 December 26 at 24 /mi, 
70/im, and 160 yum (AORKEY 4425472). Point-like sources 
were identified with the software MOPEXE] in all IRAC bands 
and in the MIPS 24 /zm band. 

The longer wavelength bands of MIPS were not used be- 
cause of the low spatial resolution. Similarly, we did not make 
use of the observation by the WISE telescope (Wright et al. 
2010) at 3.4 /zm, 4.6/im, 12/im, and 22 /mi, because of the 
coarser spatial resolution (from 6.1" to 12") as compared to 
IRAC (1.66 - 1.98") and MIPS 24/im (6"). 

The low Galactic latitude of IC 10 (-3.34°) results in a high 
extinction, E(B - V) = 0.81 (Massey & Armandroff 1995; 
Massey et al. 2007), implying that observations are strongly red- 
dened, even in the near-IR. Using the value of E(B - V) = 0.81, 
and assuming the total- to- selective extinction ratio R v = 3.05, 
we de-reddened the UBVR observations following the prescrip- 
tion of Table 3.21 in Binney & Merrifield (1998). Magnitudes 
from the / band to the IRAC [8.0] band were corrected for ex- 
tinction using the power law prescription of Martin & Whittet 
(1990), summarized in Glass (1999). 



Table 1. Candidate silicate emission sources. 



ID a (J2000), 6 (J2000) 



ID a (J2000), 6 (J2000) 



#1 

#2 
#3 
#4 
#5 
#6 
#7 
#8 
#9 



00:19:55.1, 
00:19:57.8, 
00:20:01.8, 
00:20:02.6, 
00:20:02.9, 
00:20:03.5, 
00:20:04.5, 
00:20:04.9, 
00:20:10.1, 



+59:18:34.7 
+59:18:36.6 
+59:19:34.9 
+59:17:46.6 
+59:18:27.3 
+59:18:01.4 
+59:18:50.0 
+59:18:03.3 
+59:17:39.3 



#10 
#11 
#12 
#13 
#14 
#15 
#16 
#17 
#18 



00:20:11.8, 
00:20:12.4, 
00:20:12.8, 
00:20:22.1, 
00:20:22.1, 
00:20:22.6, 
00:20:25.1, 
00:20:27.3, 
00:20:32.6, 



+59:18:25.5 
+59:17:26.3 
+59:17:10.6 
+59:17:43.9 
+59:17:24.9 
+59:17:34.6 
+59:18:08.8 
+59:15:37.0 
+59:17:05.7 



It is difficult to estimate the sample completeness based on 
the silicate strength map alone. Given the presence of gaps in 
the map (Section[2]) and given that the point spread function is 
slightly more extended than the SL slit height 3.7"), we 
could be missing somewhat less than half of the sources. The 
present sample is mostly limited by the flux at lOyum (see also 
Section|6.3|). 



3. Pixel-based analysis 

In order to identify dust enshrouded O-rich stars in IC 10, a 
silicate strength map was built from the IRS map (Section[2]). 
Silicate dust is looked for in emission via the 9.7 /im emission 
feature which originates from the Si-0 bond stretching mode 
(e.g., Knacke & Thomson 1973). Following Spoon et al. (2007), 
the silicate strength S s ii is defined as: 

c , /9.7,obs /1N 

Ssn = In- , (1) 

J9.7,cont 

where /9.7, bs is the observed flux density at 9.7 jjm, and /9.7, C ont 
is the continuum flux density at the same wavelength. We first 
calculated the continuum shape in order to estimate /9.7,cont- We 
used CUBISM to extract the continuum maps at ~ 5. 4 /mi (me- 
dian flux within 5.2-5.6/mi) and 14.0/mi (median flux within 
13.7-14.3/mi), and degraded the 5 .4 /mi map to reach the spa- 
tial resolution of the 14.0/mi map (3.7"). The continuum shape 
was then calculated by applying the spline method adapted to 
PAH-dominated spectra, as described by Spoon et al. (2007). 
Finally, we extracted the 9.7 /mi map and convolved it to a 3.7" 
resolution in order to estimate /9.7, bs- We calculated S s n using 
Equation[T] for all the pixels with a continuum flux at 9.7 /mi 
greater than 1 mJy. 

The resulting silicate emission map is shown in Figure[T] 
Based on the FWHM of sources in Figure[2| the spatial reso- 
lution ranges from 3" to 4", i.e., about 10 - 14 pc at a distance of 
700 kpc. A total of 18 point-like sources were selected based on 
their point-like appearance and for which S s n > 0.05 (Table[T} 
Figure[2]). We find no evidence of spatial clustering and no evi- 
dence of extended silicate emission. 

We wish to emphasize that the silicate strength values in- 
ferred from the pixel-based analysis (Figure[T]) are only indica- 
tive and bear systematic uncertainties. The map interpolation 
combined with the convolution to homogenize the data neces- 
sarily results in uncertain pixel values (Section|2]). Most impor- 
tantly, significant background emission arising from the ISM of 
IC 10 prevents an accurate estimate of the source silicate dust 
continuum. Final values of the silicate strength are derived using 
optimal spectral extraction (Section[5]). 

2 http://ssc. spitzer. caltech. edu/postbcd/mopex. html 



4. Optical and infrared counterparts 

4.1. 2MASS catalog 

Most candidate silicate emission sources in TableQ] could be as- 
sociated with a 2MASS source from the point-like source cata- 
log (TableQ. Given the spatial resolution of the silicate strength 
map (3" - 4"; Section|3]), we considered positive matches for 
association of 4" or less between the source centroid in the sil- 
icate strength map and the 2MASS coordinates. Most sources 
were matched within 2" or less, corresponding to about 6.8 pc in 
actual distance. There were no multiple matches within 4". We 
consider from now on the 2MASS coordinates as our reference 
coordinates for the candidate silicate emission sources. 

Sources #1, #5, #9, #17, and #18 could not be associ- 
ated with 2MASS point-like sources. Source #18 is 4.4" away 
from 2MASS 00203222+5917091 in the extended source cata- 
log. Multiple sources from the LGGS (optical) catalog are seen 
close to this s ource , but no match in the IRAC [8.0] band was 
found (Section |42| ). For this reason, and because the coordinate 
match between source #18 and the 2MASS counterpart is some- 
what larger than the map resolution 4"), we exclude this 
source from the following discussion. The match between the 
other sources and IRAC sources is discussed in Section l4~2l 

4.2. IRAC and MIPS sources 

We cross-correlated our sample with sources identified in the 
IRAC maps with MOPEX (Section[2]). At the spatial resolution 
of the IRAC [3.6] map (1.66" or * 5.6 pc), sources #4, #8, #11, 
and #12 looked by eye slightly elongated while sources #10 and 
#16 are multiple or visibly extended. Based on the radii fitted by 
MOPEX, the other sources are point-like. IRAC magnitudes are 
presented in Table[3] 

Source #17 is not a 2MASS source, but it could be matched 
with an IRAC point-like source. It is likely significantly embed- 
ded so that even the near-IR bands are extinguished by dust. 
Thus, we choose to include this source in the following discus- 
sion. Sources #1, #5, and #9 are either significantly far from any 
IRAC source or they are part of an extended emission in the 
IRAC bands. These 3 sources are excluded in the following dis- 
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35.0 30.0 25.0 20.0 15.0 10.0 05.0 0:20:00.0 19:55. 

Right ascension (J2000) 

Fig. 2. Optical image of IC 10, with I, B, and V band images as RGB colors. The outer circles (7" radius; « 24 pc) are centered on 
candidate silicate emission point-like sources (Table[T]). The inner circles (4" radius; « 14 pc) indicate the position of the 2MASS 
counterpart (Table[2]). The size of the circles is chosen for display purposes. The IRS spectra are extracted from a beam with 
FWHM- 4". 



Table 2. 2MASS Photometric data of associated stars. 



ID 


2MASS ID 


J 


H 


K 


Foreground stars 


#2 


00195768+5918349 (1.9") 


(15.19) 


(14.44 + 0.06) 


(13.71) 


#10 


00201183+5918267 (1.3") 


(13.57) 


15.39 + 0.16 


14.75 + 0.16 


#11 


00201237+5917279 (1.7") 


13.62 + 0.03 


12.56 + 0.04 


12.08 + 0.03 


#15 


00202240+5917332 (2.0") 


15.44 + 0.07 


14.60 + 0.06 


14.02 + 0.07 


#16 


00202520+5918070(1.9") 


14.95 + 0.08 


13.96 + 0.08 


(13.18) 


IC10 


#3 


00200155+5919332 (2.5") 


15.66 + 0.08 


14.60 + 0.06 


14.01 + 0.07 


#4 


00200259+5917481 (1.5") 


15.23 + 0.06 


14.15 + 0.05 


13.56 + 0.05 


#6 


00200322+5918013 (2.1") 


15.19 + 0.06 


14.05 + 0.06 


13.67 + 0.06 


#7 


00200452+5918521 (3.6") 


15.57 + 0.06 


14.4 + 0.07 


13.79 + 0.04 


#8 


00200510+5918039 (1.7") 


15.60 + 0.08 


14.81 + 0.10 


14.13 + 0.07 


#12 


00201270+5917121 (1.7") 


16.06 + 0.11 


14.88 + 0.10 


14.51+0.10 


#13 


00202225+5917432(1.4") 


15.59 + 0.09 


14.47 + 0.09 


13.76 + 0.07 


#14 


00202199+5917244 (0.9") 


14.72 + 0.05 


13.60 + 0.05 


13.24 + 0.05 



Notes. Sources #1, #5, #9, #17, and #18 could not be matched with any 2MASS sources. Magnitudes bet ween parentheses indicate upper limits 
or uncertain measurements. The field stars are identified based on the color diagnostics discussed in Section |4~4] The distance between the source 
centroid in the silicate strength map and the associated 2MASS catalog is indicated between the parentheses. 



cussion since they also do not have a 2MASS point-like source 
counterpart. 

Sources #12 and, to a lesser extent, #8 and #11 are matched 
with MIPS 24ytmi point-like sources (Figure[3]). The spatial res- 
olution at 24ytmi is 6" (corresponding to « 20.4 pc), i.e., some- 
what larger than the silicate strength map resolution (« 3.5", 
Section[3]). We verified that there were no multiple IRAC sources 
contributing to the extraction aperture at 24 yum. 



4.3. LGGS catalog 

Matches with the optical LGGS sample were found within ~ 
1" of the 2MASS coordinates for all sources in our remaining 
sampled (Table]?]) except for sources #3 and #11. Because of 
the large IRS and IRAC beams, we consider that LGGS sources 
within ~ 4" c oul d con tribute to the extracted IR fluxes. We list 
in Tables lA.ll and lA.2l the LGGS sources found within 4" of the 
2MASS counterpart. In practice, only sources #12, #13, #14, and 
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Right ascension (J2000) 

Fig. 3. Spitzer image (R: MIPS [24], G: IRAC [8.0], B: IRAC [3.6]). See Figure^for the symbol description. Sources #8, #11, and 
#12 are 24 /mi point-like sources. The gray polygon shows the area covered in all bands. 



Table 3. Spitzer photometry. 



ID 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


Foreground stars 


#2 


13.31 




13.01 






#10 a 






13.15 






#11 


11.67 


11.68 


11.35 


10.97 


8.43 + 0.30 


#15 


13.64 


13.55 


13.57 


13.15 




#16 a 


12.94 


12.87 


12.55 






IC10 


#3 


13.32 


b 


12.59 


b 




#4 


13.29 


13.16 


12.83 


12.47 




#6 


13.05 


13.04 


12.82 


12.55 




#7 


13.29 


b 


12.98 


b 




#8 


13.21 


12.65 


12.15 


11.13 


7.62 + 0.30 


#12 


13.74 


12.96 


12.23 


10.77 


5.90 + 0.10 


#13 


13.20 




12.40 






#14 


13.12 


13.16 


12.37 






Unknown membership 


#17 


14.99 


13.81 


12.90 


12.07 





Notes. Magnitudes calculated assuming zero-magnitude fluxes from the 
instruments handbooks. The uncertainties are « 0.01, « 0.01, « 0.02, 
and * 0.03 for [3.6], [4.5], [5.8], and [8.0] respectively. Bands with no 
data correspond to sources that were not found by the point source de- 
tection algorithm in MOPEX unless otherwise noted. The membership 
is based on the color diagnostics discussed in Section |4~4| 

(a) Blended or multiple objects. ^ Not covered by the IRAC 
observation. 



#16 in our sample have several bright infrared LGGS stars within 
such a radius. It must be kept in mind in the following that these 
sources could be multiple objects. 



For all the other sources, we cannot exclude that compact 
stellar clusters might be affected by confusion, even in the opti- 
cal LGGS observations. Such clusters would have to be smaller 
than ~ 1", or ~ 3pc, to be unresolved in the LGGS. We ex- 
plored the high spatial resolution observations from the Hubble 
Space Telescope (HST) to investigate further the possible pres- 
ence of compact stellar clusters. For this tes t, on ly sources that 
are member of IC 10 are considered (Section pj~4| ). Only sources 
#3, #7, #8, #12, #13, and #14 were covered by the observations 
with the Advanced Camera for Surveys (ACS) (Figure]?]). There 
is no evidence of enhanced clustering toward these sources, ex- 
cept maybe for source #8, with a few bright stars within the IRS 
extraction aperture. In the following, we assume that the flux ex- 
tracted in 2MASS, IRAC, and IRS is dominated by the LGGS 
object found closest to the 2MASS coordinates. 

4.4. Field contamination by foreground stars 

We now investigate the photometry of the objects associated 
with the candidate silicate emission sources (Tables[2]and|4]) in 
order to test their membership to IC 10 and their intrinsic stel- 
lar properties. As explained in Massey et al. (2007), the V and 
B - V colors provide a good diagnostic on the stellar type while 
also separating foreground stars from stars in IC 10. Figure|5] 
shows the photometric data from Massey et al. with the candi- 
date silicate emission sources from Table[T] overlaid. According 
to Massey et al., RSGs belonging to IC 10 are expected to have 
B-V > 2 and V < 20. Only sources #4, #6, #7, #8, #12, #13, and 
#14 fit this criterion. The other sources (#2, #10, and #15) could 
be yellow supergiants, but they are far more likely foreground 
stars. Note that sources #3 and #11 have no optical counter- 
parts (Section |43] ). Diagnostics for these 2 sources are based on 
their IR photometry alone. We also show in Figure|5]the LGGS 
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Table 4. Optical photometric data of associated stars. 



ID 


LGGS ID 


U 


B 


V 


R 


/ 


Foreground stars 


#2 


J001957.61+591835.5 (0.8") 




22.58 


20.94 


19.90 


18.77 


#10 
#11 


J00201 1.91+591827.6 (1.0") 


21.70 + 0.01 


21.29 


19.91 


19.03 


18.07 


#15 


J002022.52+591732.9 (0.9") 


19.53 


19.34 


18.40 


17.81 


17.17 


#16 


J002025.23+591807.3 (0.5") 


19.89 + 0.02 


19.96 


20.86 + 0.06 


19.45 


18.50 + 0.01 


IC 10 


#3 














#4 


J002002.61+591748.2 (0.2") 




24.67 + 0.02 


21.63 


19.62 


17.43 


#6 


J002003.23+591801.6 (0.3") 




24.02 


21.27 


19.43 


17.39 


#7 


J002004.54+591852.3 (0.2") 




25.29 + 0.08 


22.56 + 0.02 


20.35 


17.92 


#8 


J002005.11+591804.1 (0.2") 




23.35 + 0.02 


21.10 


19.43 


17.71 


#12 


J002012.73+591712.3 (0.2") 




24.54 + 0.02 


22.02 + 0.01 


20.30 


18.32 


#13 


J002022.28+591743.3 (0.2") 




23.92 + 0.01 


21.68 ±0.01 


19.81 


17.71 


#14 


J002022.01+591724.5 (0.1") 




22.74 


20.10 


18.52 


16.89 



Notes. Errors are below 0.01 dex unless otherwise noted. The field stars are identified based on the color diagnostics discussed in Section |4~4| 
The distance between the source centroid in the silicate strengt h ma p a nd the associated LGGS catalog is indicated between the parentheses. The 
LGGS ID is the best match within the search radius (see Tables [Al71 and |A.2| ). 




Fig. 4. HST/ACS images from candidate silicate emission sources within IC 10. Images were downloaded from the Hubble Legacy 
Archive (http://hla.stsci.edu/) with I, G, and B band images as RGB colors. The cross indicates the 2MASS coordinates. The circle 
represents the size of the IRS extraction aperture, 2" radius (corresponding to « 7pc). The other sources were not observed with 
ACS. 



sources found within the sea rch ra dius of the candidate silicate 
emission sources from Table |A.l| Only 3 sources have optical 
colors expected from RSGs, but they are not the most infrared 
bright within the search radius. 

In an independent approach, we have used the Besancon 
population synthesis model (Robin et al. 2003) to count the num- 
ber of foreground stars expected toward IC 10 with colors and 
magnitudes similar to those expected from RSGs in IC 10. We 
used the expected V magnitude and B-V color of RSGs in IC 10 



from Massey et al. (2007) and found that no foreground stars can 
be at the same time that red and that bright. 

Finally, near-infrared photometry further constrains the 
membership to IC 10. Based on an offset field of view, Borissova 
et al. (2000) conclude that foreground stars have colors such that 
0.4 < J - K < 1.0, and H - K < 0.1. All the sources in Table|2] 
have H - K > 0.1 and J - K between 1.4 and 1.8, which bol- 
sters our confidence that sources #4, #6, #7, #8, #12, #13, and 
#14 belong to IC 10. Moreover, Borissova et al. estimate that 
RSGs should have / - K ~ 1.4 and 13 < K < 15 while AGB 
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Fig. 5. V vs. 5 - V colors of stars toward IC 10 (Massey et al. 
2007). The candidate silicate emission sources from this study 
are shown as red points. Small black points indicate LGGS 
sources within the search radius of the candidate silicate emis- 
sion sources (Table |ATT] ). 



stars should be fainter. Source #11, which has no optical coun- 
terpart, has K = 12.1 and does not fit the constraints above; it is 
brighter in the K band by one order of magnitude than the other 
sources and is likely a foreground star. All the other sources have 
1.4 < J - K < 1.8. In particular, source #3, which also has no 
optical counterpart, could be a RSG in IC 10 based solely on its 
JHK colors. 

In summary, both the optical and IR photometry of sources 
#4, #6, #7, #8, #12, #13, and #14 are consistent with member- 
ship in IC 10. Furthermore, the K magnitude of these sources 
(and of source #3) seem to indicate they are RSGs. We refine the 
determination of the stellar nature in Section[6] 



5. Mid-infrared spectra 

In this section, we present the mid-IR Spitzer/IRS spectra of all 
the candidate silicate emission sources. We consider sources #3, 
#4, #6, #7, #8, #12, #13, and #14 (members of IC 10) as well as 
#17 (IRAC source). 

The presence of spatially extended MIR emission (domi- 
nated by polycyclic aromatic hydrocarbon bands and warm dust 
continuum) from the ISM of IC 10 prevents a regular spectral 
extraction of the sources (i.e., integration of the flux within a 
spatial window). We therefore used the optimal extraction pro- 
vided by SMART- AdOpf] (Lebouteiller et al. 2010) to extract 
the spectra at the matching star location in the exposure images. 
Optimal extraction weighs the source spatial profile by using the 
instrument point spread function as a reference. The extended 
emission coming from the ISM of IC 10 was removed simul- 
taneously using a second- or third-order polynomial. The stars 
location within the IRS SL aperture was constrained not only in 
the cross-dispersion direction but also in the dispersion direc- 
tion, therefore accounting for the slit throughput and providing 
an accurate flux calibration. An example of extraction is pre- 
sented in Figure|6] Final spectra are presented in Figure[7] 

The spectral trace was detected for all the sources mentioned 
above. Detection levels (based on the integrated SL wavelength 
range) are given in Table|5] Besides the bright sources #8 and 



#4- 




8(jm 



11.3|il71 



14|am 




w (cross-dispersion coordinates in detector pixels) 

Fig. 6. Example of optimal spectral extraction (source #4). The 
detector image is shown in the top panel, after background sub- 
traction, with the cross-dispersion profile (w) as a function of 
wavelength. The corresponding profile along the aperture (inte- 
grated signal over the wavelength range vs. w) is shown as a 
histogram in the bottom panel. The connected squares show the 
fit of the spatial components in the slit, including source #4 (red 
profile), the extended background emission, and another slightly 
extended source in IC 10 matching the location of the H n re- 
gion [HL90] 17 (Hodge & Lee 1990), here fitted with 2 point- 
like sources showed by the green and blue profiles). 



#12, we note that sources #4, #6, #13, and #17 are fairly well 
detected (more than 2<x), while sources #3, #7, and #14 barely 
stand above the detection threshold. Based on the comparison 
between the source spatial profile and the IRS point spread func- 
tion, we find that all sources are point-like at the spatial resolu- 
tion of the IRS SL module at 10/mi, i.e., « 2". 

Sources #8 and #12 show prominent silicate emission peak- 
ing respectively at 8 mJy, and 1 1 mJy, while source #4 shows 
weak emission peaking at « 3mJy (Figure[7]). The signal-to- 
noise ratio of the other sources is too low to assert unambigu- 
ously the presence of silicate dust. Silicate strength values are 
given in Table[5] The silicate emission flux density peak in RSGs 
of the Large Magellanic Cloud (LMC) ranges from ~ 1 Jy to 
~ 3.5 Jy (Buchanan et al. 2009). Such sources would have fluxes 
around 4 - 14mJy at the distance of IC 10, which is compati- 
ble with our values. We notice that source #12 is characterized 
by a bright dust continuum longward of 13 jim which cannot be 
due to background extended emission, as it was removed dur- 
ing spectral extraction. This object is associated with a point- 
like source in the MIPS 24 yum image. The nature of the silicate 
emission sources is discussed in SectionETTl 



Version 8.2.4; http://isc.astro.cornell.edu/IRS/SmartRelease 



6. Properties of the stars 

6.1. RSG vs. AGS 

Silicate dust can be produced by several types of stars, most no- 
tably AGBs, RSGs, planetary nebulae, and novae. Young stellar 
objects (YSOs) also show silicate dust in their disks or outflows, 
although the dust grains might not be produced in situ. The opti- 
cal and near-IR colors suggest that all sources but #17 could be 
RSGs (Section 4.4). We now review this finding by comparing 
the K s ws. J-K s colors of stars in IC 10 and in the LMC. Figure[8] 
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6 8 10 12 14 



Wavelength [pm] 

Fig. 7. Spitzer/IRS spectra of sources photometrically identified as IC 10 sources. The 3 filled circles represent the observed con- 
tinuum flux density at 8.0 /zm, 9.7 /mi, and 14.0 /mi while the open circle represents the interpolated flux density at 9.7 yum used to 
infer the silicate strength (see text). Spectra were smoothed by a running 3 -pixel median box. 



Table 5. Silicate strength values and model results. 



ID 


Detection 


S sil 


^bol 


?eff 


MLR^ 






(<r) 


(mags) 


(L ) 


(K) 


(Moyr" 1 


) 


#3 


1.4 


+ 1.42 +a63 
+2.22+j : § 


120000 


3 550 


8 x IO" 7 




#4 


3.4 


130000 


3 397 


2 x IO" 7 




#6 


2.5 


-0 53 +0 ' 75 

- a89 S 

+ 1 44+0:07 


130000 


3 490 


2 x IO" 7 




#7 


1.5 


110000 


3 397 


3 x IO" 7 




#8 


7.7 


120000 


3 550 


28 x 10" 


7 


#12 


18 


+0.98^ 


90000 


3 550 


30 x 10" 


7 


#13 


2.4 




110000 


3 550 


14 x 10" 


7 


#14 


1.3 


TZ,.JJ_2 49 


150000 


3 550 


5 x IO" 7 




#17 


2.3 


+ 1 27 +o:85 


15 000 


3 550 


30 x 10" 


6 



Notes. The sources with significant silicate emission are shown in bold 
(Section[5j. 

(a) Detection level over the spectral trace (integrated SL wavelength 
range). {b) Mass-loss rate determinations with a factor of « 2 statistical 
uncertainties. 



shows the CMD in which the magnitudes of IC 10 sources have 
been scaled to the distance of the LMC. Sources #7 and #13 ap- 



pear to fall in the O-rich AGB color domain, while source #3 lies 
on the AGB/RSG cut. All the other sources (#2, #4, #6, #8, #12, 
and #14) are unlikely to be AGB stars. 



The most important constraint in distinguishing between 
AGBs and RSGs is the bolometric luminosity. It is ex- 
pected that RSGs have a bolometric luminosity Mb i < -7.9 
(i.e., >117 000L o ), while AGBs should have M bol > -7.1 
(<56000L o ) (e.g., Wood et al. 1983). Although AGBs can un- 
dergo hot bottom burning or thermal pulses that can increase 
their brightness temporarily (e.g., Groenewegen et al. 2009), the 
threshold M\> \ < -7.9 allows separating less luminous RSGs 
from intermediate-mass AGBs (see Massey et al. 2003; Massey 
& Olsen 2003). We used the MIR spectra alongside the photom- 
etry to constrain the luminosity and mass-loss rate of the sources 
using the radiative transfer model described by Groenewegen et 
al. (1995, 2009). For all stars we fitted a model with pure silicate 
dust (with absorption coefficients from Volk & Kwok 1988), and 
another one with a mixture of 20% aluminum oxide and 80% 
silicate. For #8 and #12 the pure silicate dust model provided 
the best fit, while either model fits the data for the other sources. 
Tableland Figure[9] shows the results for the silicate dust model. 
Based on the model results, most of the sources are much too 
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J-Ks 

Fig. 8. ^ s vs. J- K s CMD for sources in IC 10. The background 
is a Hess diagram of the sources from the SAGE catalog of the 
LMC (Meixner et al. 2006). The IC 10 sources are corrected 
for reddening and then "moved" to the LMC for comparison. 
AGB stars fall within the labeled region in the upper-right region 
(Cioni et al. 2006; Nikolaev & Weinberg 2000). O-rich AGB 
stars are redder than the left oblique line while C-rich AGB stars 
are redder than the right oblique line. Sources #10, #11, #15, 
and #16 are foreground sources not associated with IC10, and 
are marked with green triangles. Source #17 is not a 2MASS 
source. 



bright to be AGBs. Only sources #12 and #17 fall below the lu- 
minosity threshold. Sources #3, #4, #6, #7, #8, #13, and #14 are 
thus again compatible with RSGs. 

The specific case of source #12 is puzzling, as it shows 
strong silicate emission, a dust continuum longwards of lOyum 
which is visible in the IRS spectrum, and it is detected in 
the MIPS 24 /mi image. The Herschel/PACS observations (S. 
Madden, private communication) place strict upper limits on the 
dust emission at far-infrared wavelengths. Due to the presence 
of a dust continuum, we compared the spectrum of this source to 
young star cluster dust radiative transfer models (Whelan et al. 
2011) in order to test the hypothesis that source #12 is a young 
compact stellar cluster. If source #12 is a young stellar clus- 
ter, then the near-IR measurements could arise from the main- 
sequence stars in the cluster while the long- wavelength contin- 
uum would arise from dust heated in the intra-cluster medium 
by the young stars. However, no models fit the data points: the 
24 fim data point and the PACS upper limits were much lower 
than expected for a range of appropriate dust geometries, from 
optically thin and geometrically thick to optically thin and geo- 
metrically thin. This suggests the lack of a carbonaceous grain 
dust component as one would expect in a dusty star forming en- 
vironment. This finding is compatible with the lack of a cluster- 
ing around source #12 in the HST images (Figure)?]). We there- 
fore tentatively conclude that source #12 is not a young compact 
cluster but is a dust-producing evolved star of some kind. Due to 
the presence of a silicate emission feature, low dust continuum, 
and relatively low luminosity (compared to RSGs) of planetary 
nebula NGC 6804, source #12 may be a planetary nebula around 
an O-star (Bilikova et al. 2012; Weidmann & Gamen 2011), or 
a dusty WR star. The shallow long-wavelength SED is also rem- 
iniscent of extreme- AGB stars, which are usually carbon stars 
(see Boyer et al. 2012 and references therein). 



6.2. Mass-loss rates 

Although it is possible to infer the mass-loss rate from molecular 
infrared transitions (Matsuura et al. 2006), the dust emission as- 
sociated with the circumstellar envelope is the best tracer of mass 
loss. The mass-loss efficiency in O-rich stars depends more on 
metallicity than in C-rich stars, because O-rich dust depends on 
metallicity-limited elements (Si, Al), while amorphous carbon 
depends on self-produced C (e.g., Sloan et al. 2008, 2012). 

The dust production rate was computed from our models 
(Section 



6.1 



assuming an outflow velocity of lOkms -1 . The 
total mass-loss rate is estimated using a standard gas-to-dust ra- 
tio of 200. The mass-loss rates we infer (Table[5]) lie within the 
ranges of what is found for the RSGs in the SMC and LMC, with 
rates between 10" 5 5 M© yr" 1 and 10" 7 M© yr" 1 (Groenewegen et 
al. 2009, assuming identical values for the outflow velocities and 
gas-to-dust ratio). 



6.3. Spatial distribution and sample completeness 

A comparison of the luminosities of the RSG candidates 
(Table|5]) and theoretical isochrones (Fagotto et al. 1994) sug- 
gests that the stars in our sample (all but #17) are more massive 
than > 12 M Q . Their expected lifetime is ~ 20 Myr old, which is 
consistent with a starburst population. Is the spatial distribution 
of the RSG candidates compatible with the starbursting region? 
The claim for a starburst in IC 10 mainly originates from the dis- 
covery of over 100 WR stars (Massey et al. 1992; Royer et al. 
2001; Massey & Holmes 2002). These studies showed that the 
spatial distribution of the WR stars is quite uniform, suggesting a 
widespread starburst. Hence, we do not expect to find the RSGs 
in any particular region, which is supported by our results. 

We expect the number of RSGs to dominate over the number 
of WR stars at low metallicity, with a lower limit of 50 RSGs 
in IC 10 (e.g., Maeder et al. 1980; Massey et al. 2002, 2003). 
Our current sample sets a lower limit on the actual number of 
O-rich dust enshrouded stars and RSG candidates, with notably 
a strong limitation by the Spitzer/IRS sensitivity at lOyum. In 
order to quantify the completeness of the sample, we built a 
sample of stars with similar colors as the RSG candidates we 
already identified. Based on Sections |4~4| and |6.1| we choose 
the following constraints: J - K = 1.5 ± 0.2, H - K > 0.3, 
B - V > 2.2, and V - R > 1.6. Only 16 sources in the 
2MASS/IRAC cross-matched sample (633 stars) fit these con- 
straints, including the already confirmed sources #4, #8, and #12 
(Table|6]). None of the other 13 sources show silicate emission 
in their IRS spectra. We partly attribute this low number of sili- 
cate emission sources to the low signal-to-noise ratio in the IRS 
spectra, as indicated by the IRAC [8.0] magnitudes. Only one 
source, 2MASS 00200459+5918198, is expected to be bright 
enough for the silicate emission to be detected, and its [5. 8] -[8.0] 
color does suggest the possible presence of silicate emission. 
However, the [8.0] flux is likely overestimated due to a contami- 
nation by PAH emission which might not have been completely 
subtracted when performing the aperture photometry with a sky 
annulus. A total of 13 more sources could thus also be RSGs in 
addition to the 3 sources we already identified. Co nsid ering an 
average mass-loss rate of 5 x 10" 7 M yr -1 (Section 6.2 ), this re- 
sults into a total mass -los s rate of 8 x 10" 6 M© yr" 1 for IC 10 (see 
discussion in Section |6~4| ). 

The low number of RSGs was already noticed by Massey 
et al. (2007), using deep optical images (photometry errors of 
0.004 in V and 0.015 in B for B = 24.3) and better spatial reso- 
lution than 2MASS and Spitzer observations. The authors argue 
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Fig. 9. SED fit results. For each source, the SED model is shown on top as the solid curve, with the diamonds indicating photometry 
points from the optical bands, 2MASS, and Spitzer, and with the segments indicating the IRS data. The [24] data point for source 
#8 bears significant uncertainties (see Table[5]). The downwa rd ar rows for source #12 show the upper limits on the Herschel/PACS 
fluxes. Source #12 could not be fitted by our models (Section[6J]). 



that a very young burst (< lOMyr) could be responsible for the 
large WR/RSG population ratio. Although our Spitzer data un- 
cover just the tip of the RSG iceberg, the missing RSGs in IC 10 
remains an unsolved mystery. 



6.4. Discussion 

Our analysis of stellar properties shows that the luminosities of 
red-supergiants in IC 10 are above 90, 000 L (Table 7). Stellar 
evolution models show that stars with an initial mass higher than 
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Table 6. Sources with similar colors as confirmed silicate emis- 
sion sources. 



2MASS source 


K 


[8.0] 


[5.8H8.0] 


00201538+5919070 


13.52 






00200259+5917481 (#4) 


13.57 


12.47 


0.36 


00200322+5918013 


13.67 


12.55 


0.27 


00202465+5919003 


13.92 






00203019+5917154 


14 






00200459+5918198 


14.09 


12.45 


0.78 


00200510+5918039 (#8) 


14.13 


11.13 


1.02 


00200825+5919092 


14.26 






00200277+5917564 


14.28 


13.21 


0.21 


00202179+5917477 


14.36 






00195375+5918118 


14.43 






00202036+5918205 


14.51 






00200839+5916419 


14.51 


14.07 


0.41 


00201270+5917121 (#12) 


14.51 


10.77 


1.46 


00200848+5916552 


14.76 


13.65 


0.37 


00200819+5919202 


15.06 







Notes. Constraints on the colors are / - K = 1.5 + 0.2, H - K > 0.3, 
B-V >2.2, mdV-R> 1.6. 



11.7M© can reach luminosities higher than 9O,OOOL during 
the RSG phase (Fagotto et al. 1994). Stars lower than 9 M Q mass 
cannot reach such a high luminosity, though the calculated mass 
range lacks 9 - 1 1.7 M stars, which evolve into the super- AGB 
phase. The models further predict that the age to reach such high 
luminosity phase is about 20 Myrs old or younger. 

The age of RSGs is consistent with the star-formation his- 
tory of this galaxy. Hunter (2001) analyzed stellar clusters in 
the galaxy, and showed that IC 10 has episodes of high star- 
formation; young (4 - 30 Myrs) clusters, presumably formed in 
the starbursts, and intermediate age (450 Myrs) clusters. The age 
of RSGs corresponds to the starburst phase of this galaxy and 
it is consistent with the large population of WR stars detected 
(Massey & Holmes 2002). 

Our analysis might provide a unique case of measuring 
mass-loss rates of such young red-supergiants. RSGs in IC 10 
might belong to young starbursts (4 - 30 Myrs), while LMC 
RSGs are mostly from intermediate age clusters (a few 100s of 
Myr; Elson & Fall 1988; van Loon et al. 1999). The mass-loss 
rates of young RSGs are typically 10" 7 - 10" 6 M© yr" 1 , which are 
comparable to those observed in LMC RSGs. We did not detect 
mass-loss rates higher than 10" 6 M yr -1 in our sample of RSGs 
(if removing #17), although such higher mass-loss rates have 
been found in Galactic RSGs (e.g., VY CMa of 10" 4 M Q y^; 
Decin et al. 2006). This might be due to the fact that our se- 
lection of RSGs is limited by the detection limit of the 2MASS 
photometry, where RSGs with high mass-loss rates are faint at 
near-IR wavelengths. 

The current analysis of mass-loss rates shows that gas ejected 
from the RSG population into the ISM is at least 10" 5 M© yr" 1 , 
and the dust return from RSGs is at least 5 x 10" 8 M© yr" 1 . We 
consider this to be a lower limit, because it is likely that stars 
with high mass-loss rates are not detected yet (because of dust 
extinction in the K band), but these stars could to contribute 
significant fraction to gas and dust return from the total RSG 
populations into the ISM. Compared with the gas mass of IC 10 
(~ 10 8 M Q , Yin et al. 2010) and a star-formation rate of up to 
0.2M o yr _1 (Leroy et al. 2006), which represents the ISM gas 
mass consumed by the formation of stars, the mass injected from 
RSGs is significantly small, and has little impact on the total gas 



ISM mass at the current stage. Although the total mass of ISM 
dust in this galaxy is unknown, it is likely to be on the order of 
1O 6 M , considering the gas-to-dust mass ratio. The dust from 
RSGs appear to be not an important contributor to the dust mass 
in this galaxy if the star formation rate of this galaxy has been 
more or less similar for the past few Myrs. A similar conclusion 
was reached for the Magellanic Clouds where AGBs are shown 
to dominate the stellar dust production (e.g., Boyer et al. 2012). 
We conclude that although RSGs could in principle dominate the 
dust production over AGBs in a starbursting galaxy (e.g., Massey 
et al. 2005), this is not observed in IC 10. 

7. Conclusions 

We report the discovery of O-rich dust enshrouded stars within 
the nearby 700 kpc) dwarf starburst galaxy IC 10. We exam- 
ined the Spitzer/IRS spectral map (7.5 - 14.5 yum) in order to 
build a sample of point-like sources potentially showing silicate 
dust in emission. The silicate strength map we constructed re- 
veals several point-like sources and no extended emission. 

Most sources are associated with single, point-like, 2MASS 
and optical sources. We investigate the colors and magnitudes in 
the near-IR and optical, and identify 9 sources belonging to the 
IC 10 system. The colors and photometry in the optical and near- 
infrared suggests that these sources are different from AGB stars. 
Modeling of the dust predicts high luminosities compatible with 
RSGs. This is thus the farthest detection of O-rich dust encircled 
stars confirmed spectroscopically. The low number of sources 
discovered spectroscopically does not solve the apparent lack of 
RSGs as compared to WR stars in IC 10 (Massey et al. 2007). 

We derived mass-loss rates for all sources using a radiative 
transfer model. Accounting for sample completeness, the total 
mass-loss rate is significantly small as compared to the dust mass 
in IC 10. Other sources of dust (AGBs, SNe, WR) are necessary 
to explain the dust mass observed in the ISM of IC 10. 

Another source belonging to the IC 10 system shows strong 
silicate emission together with a warm carbonaceous dust grain 
continuum. The nature of this source remains unknown. 
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Table A.l. LGGS sources within d < 4" of the silicate emission 
candidates identified as members of IC 10. 



ID 


d(") 


LGGS 


V 


V - R 


R- I 


#3 


2.33 


J002001. 84+591933.9 


23.463 


0.684 






2.89 


J002001.22+591934.6 


22.928 


0.647 


1.032 




3.40 


J002001. 82+591930.5 


22.891 


0.691 


0.653 




3.83 


J002001.96+591935.4 


23.475 


1.025 


0.904 


#4 


0.18 


J002002.61+591748.2 


21.629 


2.008 


2.192 




2.22 


J002002.55 +591750.3 


23.484 


1.020 


1.241 




3.29 


J002002.94+591746.2 


21.804 


0.492 


0.480 




3.79 


J002003.05+591749.5 


22.522 


0.440 




#6 


0.31 


J002003.23+591801.6 


21.274 


1.848 


2.035 




2.27 


J002003 .48 +59 1 800.2 


23.084 


0.537 






2.59 


J002003.31+591803.8 


23.142 


0.938 


1.183 




3.77 


J002002.94+59 1804.4 


21.656 


0.745 


0.848 




3.91 


J002003.25+591805.2 


23.494 


0.703 




#7 


0.25 


J002004.54+591852.3 


22.559 


2.207 


2.428 


#8 


0.21 


J002005.11+591804.1 


21.096 


1.666 


1.717 




0.98 


J002005. 19+591804.6 


21.529 


0.238 


1.206 




2.52 


J002005.00+591806.3 


22.979 


0.720 


0.601 




3.60 


J002005.45+591801.5 


19.672 


0.875 


0.997 


#12 


0.30 


J002012.73+591712.3 


22.020 


1.722 


1.978 




2.23 


J002012.75+591709.9 


22.271 


1.535 


1.453 




2.51 


J002013.00+591713.1 


23.496 


0.845 


1.139 




3.02 


J002012.59+591715.0 


22.626 


0.689 


1.024 




3.64 


J002012.90+591715.4 


23.185 


0.928 


1.220 




3.82 


J002013. 19+591711.4 


21.644 


0.909 


0.976 


#13 


0.25 


J002022.28+591743.3 


21.682 


1.870 


2.095 




1.54 


J002022.45+591743.0 


22.817 


0.632 






2.49 


J002022.46+591741.3 


21.704 


0.358 


0.290 




2.93 


J002021.97+591741.2 


20.790 


0.492 


0.141 




3.30 


J002022.23+591746.5 


22.346 


2.662 


1.090 




3.33 


J002022. 13+591746.4 


20.599 


1.013 


0.972 




3.39 


J002022.68+591742.4 


21.850 


0.478 


0.779 




3.59 


J002021. 87+591741.1 


19.094 


0.452 


0.447 


#14 


0.18 


J002022.01+591724.5 


20.096 


1.581 


1.626 




0.84 


J002022.05+591725.1 


22.266 


3.631 


1.774 




1.63 


J002022.10+591725.8 


21.842 


1.368 


0.601 




1.76 


J002022.1 1+591725.9 


22.051 


1.566 


0.616 




1.89 


J002021.78+591725.4 


22.936 


0.132 






2.25 


J002022.22+591723.0 


23.046 


1.235 


1.264 




2.35 


J002022.05+591722.1 


21.650 


0.932 


1.078 




2.45 


J002021.85+591722.2 


22.018 


0.846 


1.099 




2.82 


J002021.82+591726.9 


23.105 


1.172 


1.418 




3.75 


J002021.68+591721.5 


20.448 


0.514 


0.444 




3.94 


J002022.44+591726.3 


22.516 


0.469 


0.702 




3.96 


J002021.60+591721.8 


20.446 


0.512 


0.444 



Notes. The LGGS ID in bold indicates the best match in terms of dis- 
tance and infrared brightness. 



Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 
Yin, J., Magrini, L., Matteucci, E, et al. 2010, A&A, 520, A55 

Appendix A: Cross-correlation LGGS stars - 
silicate emission candidates 

The silicate emission candidates from TableQ]are matched with 
optical sources from the LGGS in Section |4~3| When several 
LGGS stars fall within 4" of the silicate emission source coor- 
dinates, we selected the closest LGGS star. We list in Tables|A3] 
and IA.1I the LGGS sources within 4" from each silicate emis- 
sion candidate. We also list their optical colors. In all cases, the 
closest LGGs star is also the best match in terms of brightness 
and red color. 
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Table A.2. LGGS sources within d < 4" of the silicate emission 
candidates identified as field stars. 



ID 


</(") 


LGGS 


V 


V-R 


R-I 


#2 


0.80 


J001957.61+591835.5 


20.938 


1.036 


1.132 




3.12 


J001957. 83+591 837.8 


23.890 


0.997 


1.084 




3.49 


J001957.24+591835.8 


22.182 


1.589 


1.717 




3.59 


J001957.99+591832.2 


22.946 


0.550 


0.604 


#10 


1.09 


J002011.91+591827.6 


19.911 


0.883 


0.955 




2.51 


J002011.55+591825.4 


21.903 


0.672 






3.02 


J002011.87+591829.7 


22.798 


0.438 






3.74 


J00201 1.76+591823.0 


22.901 


1.456 


1.408 




3.95 


J002011.59+591823.2 


23.305 


1.091 


1.452 


#11 


3.78 


J002012.61+591724.6 


23.381 


0.834 


0.997 


#15 


0.97 


J002022.52+591732.9 


18.404 


0.590 


0.642 




1.81 


J002022.17+591732.8 


22.633 


1.026 


1.019 




2.46 


J002022.33+591730.8 


22.744 


1.521 


1.460 




3.52 


J002021.99+591731.6 


23.073 


0.443 


0.921 




3.60 


J002022.40+591729.6 


23.417 


1.003 


0.782 




3.82 


J002022.20+591729.7 


23.299 


-0.084 




#16 


0.38 


J002025.23+591807.3 


20.857 


1.409 


0.949 




0.46 


J002025.14+591807.0 


19.644 


0.263 


0.412 




1.29 


J002025.36+591806.6 


20.557 


1.082 


1.644 




1.40 


J002025 .05 +59 1 807 . 8 


20.807 


0.459 


0.188 




1.64 


J002025.41+591806.7 


21.316 


0.165 






1.92 


J002025.24+591805.1 


19.392 


0.712 


0.724 




1.97 


J002025.33+591808.7 


21.954 


0.427 






2.04 


J002024.95+591806.3 


19.696 


0.437 


0.408 




3.14 


J002025.60+591807.7 


21.873 


0.433 


0.415 




3.39 


J002025.45+591804.2 


21.704 


1.354 


1.352 




3.43 


J002025.04+591810.2 


22.570 


-0.150 






3.48 


J002025.47+591804.2 


21.753 


1.400 


1.355 



Notes. The LGGS ID in bold indicates the best match in terms of dis- 
tance and infrared brightness. 
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